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Abstract. Recently, there has been suggestions that the apparent accelerated ex- 
pansion of the universe is due not to a cosmological constant, but rather to inhomo- 
geneities in the distribution of matter. In this work, we investigate a specific class 
of inhomogeneous models that can be solved analytically, namely the dust-dominated 
Lemaitre-Tolman-Bondi universe models. We show that they do not permit acceler- 
ated cosmic expansion. 
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1. Introduction 

Kolb and coworkers pQ argued recently that growing perturbations on a scale larger 
than the Hubble length may lead to accelerated expansion as observed from the center 
of the perturbation. This was then followed up by Wiltshire [2j and Carter et al. [Sj, 
who considered a dust dominated universe model where the observed universe is an 
underdense bubble in an Einstein-de Sitter universe. They calculated the luminosity- 
redshift relationship for their model and found very small deviations from the standard 
ACDM model. Thus, they were able to obtain a good fit to the Hubble diagram of type 
la super novae. 

However, critical papers have appeared. We will now proceed to give a brief 
overview over some of these papers and the arguments therein. By considering a 
universe model where sub-Hubble perturbations are absent, Flanagan [3] argues that 
the contributions from the super-Hubble perturbations to the value of the deceleration 
parameter are so small that they cannot be responsible for the acceleration of the 
universe. Geshnizjani et al. |2j argue that to second order in spatial gradients the 
super-Hubble perturbations only amount to a renormalization of local spatial curvature, 
and thus cannot account for the negative deceleration parameter. Making an expansion 
to Newtonian order in potential and velocity, but taking into account fully non-linear 
density inhomogeneities, Siegel and Fry [0] obtain similar results, and conclude that 
inhomogeneity contributions cannot mimic the effects of dark energy or induce an 
accelerated expansion. 

One of the strongest and most convincing critiques of the conclusions of Kolb et al. 
is that of Hirata and Seljak jjj. They argue from the Raychaudhuri equation that in a 
dust dominated universe there must be a non-vanishing vorticity in order to obtain a 
negative deceleration parameter. Then they showed that the perturbations considered 
in pp have vanishing vorticity and hence cannot lead to accelerated expansion. In yet 
another critique of the work of Kolb et al., Rasanen [H] also claims that their conclusions 
are ruled out. The basis for this conclusion was an analysis of the model proposed in 
by applying the Buchert formalism for backreaction 9J. This conclusion is also 
supported in works by Giovannini jTOJE], where the author performs a general analysis 
of the class of models considered in pQ. He finds that dust-dominated models cannot 
have accelerated expansion, and hence, the conclusion of Kolb et al. cannot be correct. 

Several of the arguments against accelerated expansion induced by inhomogeneities 
in matter are based upon perturbation calculations. However, as noted in one way 
of evading these arguments is through non-perturbative effects. This is what we will 
investigate in the present work. We want to know if Einstein's field equations permit 
inhomogeneities to change the positive deceleration parameter of a dust-dominated 
homogeneous universe model to a negative value. 

Unfortunately, a general analysis incorporating all possible inhomogeneous models 
is not feasible, so we will instead concentrate on one specific class of such models, 
namely the dust-dominated, spherically symmetric models. Moffat suggested in 



Can a dust dominated universe have accelerated expansion? 



3 



that it might be possible to find solutions with accelerated expansion. Considering the 
field equations for these models, an expression for the local deceleration parameter was 
derived, and it was claimed that this might indeed allow for negative values. Later, 
in this claim was moderated considerably. Here it was claimed that the volume 
averaged deceleration parameter, and not the local deceleration parameter, could be 
negative. No explicit solution of the field equation were actually found in these analyses. 
Instead, the conclusions were reached by investigating the field equations directly. 

In this work, we will consider the same class of spherically symmetric, 
inhomogeneous universe models and present a general solution to the field equations. 
We then show that no such model can have an accelerated expansion in the sense of a 
negative local deceleration parameter. 

The structure of this paper is as follows. In Sect.|2]we present the field equations for 
a general spherically symmetric, inhomogeneous universe model. In Sect. El we discuss 
the properties of these models and whether they can have accelerated expansion. Finally, 
in Sect. E]we give a brief summary of our work and present our conclusion. 



2. Spherically symmetric, inhomogeneous universe models 

The spherically symmetric, inhomogeneous universe models are described by the 
Lemaitre-Tolman-Bondi (LTB) space-time E3- The line element can be written 

as 

ds 2 = -dt 2 + X 2 {r, t)dr 2 + R 2 {r, t)dQ 2 , (1) 

where X(r,t) and R(r,t) are general function to be determined by the field equations 
and boundary and initial conditions. They can be thought of as generalized, position 
dependent scale factors in the radial and transverse directions. The Einstein field 
equations are 

G[iv = R^iv ~ ~^9iivR = kT^v , (2) 
where k = 8nG. 

We assume the universe to be occupied by an ideal fluid with energy density p(r, t) 
and an isotropic pressure p(r,t). In the comoving coordinates defined in Eq. (fl|). the 
energy-momentum tensor of the fluid can be written as T^ u = diag (p,p,p,p). A relation 
between X and R can be found by solving the t-r component of the field equations, for 
which the right-hand side vanishes. The relation one finds is 

*M) = -^p (3) 

where f(r) is an arbitrary function of r only. Throughout this paper, we will use a 
' = d/dr to denote differentiation with respect to r and ' = d/dt for differentiation with 
respect to t. 

Following Moffat O HB] , we define two " Hubble parameters" 

tf± = f and tf r = §, (4) 
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which are a measure of the local expansion rate in the perpendicular and radial directions 
respectively, and the deceleration parameter 

Using these definitions, the time-time and space-space components of the field equations 
take the form 

and 

- 6H 2 ± q ± + 2H 2 ± - 2 _ 2H r H ± + = ~<P + 3p) • (7) 

where (3 = f 2 — 1. These are the generalization of the Friedmann equations in a 
homogeneous universe to a spherically symmetric, inhomogeneous universe. Note that 
the function (3{r) is often written as j3{r) = 2E(r)r 2 , where E(r) determines the local 
curvature radius. 

3. Expansion rate of dust-dominated LTB models 

From now on we will assume the cosmic fluid to be pressure-less matter. Adding Eqs. © 
and (0), we then obtain the following expression for the deceleration parameter 

The condition for accelerated expansion now takes the form 

P > R 2 > or f 2 > 1 + R 2 . (9) 

Substituting the definition for the deceleration parameter in Eq. (jSJ) into Eq. (JHJ), we 
arrive at the expression 

2RR + R 2 = /3. (10) 

A single integration of this equation yields 

RR 2 = PR + a or H 2 ± = A + ^ , (11) 

where the function a(r) enters as an integration "constant" and depends only on the 
radial coordinate r. Comparing this equation with the ordinary Friedmann equation for 
homogeneous models, we see that the dynamical effects of j3 and a are similar to those 
of curvature and dust, respectively. Therefore, a(r) is regarded as the gravitational 
mass function, and one often chooses a(r) cx r 3 . Note, however, that a does not appear 
in the energy-momentum tensor. 

From Eqs. © and (fTTj) . we find that a and (3 must satisfy the inequality 

-{3R<a(r)<0 (12) 

in order for the deceleration parameter to be negative. Thus, it appears that the LTB 
models seem to allow for accelerated expansion even for dust-dominated universe models, 
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as long as a and (3 are chosen appropriately. The dynamical effect of a would then 
correspond to that of dust with negative energy density in a homogeneous universe 
model. It should be noted that the inequality above forbids accelerated expansion in a 
"big bang" model where the scale factor has the initial value R(0, r) = 0, which implies 
that a(r) > 0. However, this initial condition may not be physically realistic. The 
universe may have started with a finite scale factor, or maybe has collapsed and reached 
a finite minimum radius before expanding again. In such models accelerated expansion 
does not seem to be forbidden. However, as we will show shortly, such solutions are 
unphysical. First, we find a general solution of the field equation and show how such 
an accelerated solution behaves. 

We introduce a conformal time 77 defined via the differential relation (3 l l 2 dt = Rdr]. 
This allow us to integrate Eq. (|TT| to give a parametric solution in terms of rj. The 
solutions can be grouped into three different classes according to the sign of the local 
curvature (3. Choosing t = 7/ = 0, the three classes are 



R 



a: 



(cosh?7 — 1) + R 



(3t 



^(sinliT/ 



for positive (3, 



R 



a 



2|/3| 



a 



for negative (3, and finally 



'1 - COST]) + R 



(rj - sin n) + R 



cosh T] + 



sinh rj + 



cos 77 + 



sin 77 + 



'a + (3Rq 



PRo 



sinh?7 



la + pR 
(3Ro 



(cosh 77 — 1) 



(13) 

(3 > 



a + PR 
I \P\Ro 
la + f3R c 
W\Rn 



sin 77 



^1 — cos 77) 



(14) 



< 



R 



2/3 



(15) 



for vanishing (3. We allow the initial size of the universe to be non-zero, hence the non- 
vanishing value for R = R(t = 0,r). Eqs. (|T3|l - (fT5|l with i? = represent the usual 
form of the LTB solution of the field equations. The scale factor R as a function of time 
is shown for some typical parameters in Fig. ^ The top curve represents a solution with 
accelerating expansion. 

In order for a dust-dominated solution in the LTB space-time to undergo an epoch 
of accelerated expansion, the inequality in Eq. ()12j) has to be satisfied. Can a physically 
realistic solution satisfy this? Differentiating Eq. (fTTj) and comparing the result with 
Eq. © we obtain a relation between the density distribution and a for this class of 
models, 
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where V = R 3 is a comoving volume. A physically realistic model must have positive 
density everywhere, and hence, it must have a' > 0. Furthermore, the angular part of 
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Figure 1. The scale factor for some typical values of a and (5 



the line element is R 2 dQ 2 where dVL is a solid angle element. At the origin, r = 0, this 
part of the line element must vanish, and thus R(0,t) = 0. From Eq. we then get 
that a(0) = 0. Since a'(r) > 0, it follows that a(r) > for all r. But, as we showed 
above, accelerated expansion is only possible for models with a{r) < 0. Hence, the dust 
dominated LTB universe models must have decelerated expansion. 

We shall now arrive at the same conclusion by an alternative approach related to the 
Raychaudhuri equation, which describes the flow of a cosmological fluid in space-time. 
It can be written as 



f)2 

+ j + 2(o- 2 



(17) 

is the expansion rate of 
the fluid and a and to are the shear and vorticity scalars, respectively. The latter two 
quantities are defined as 



where a A 



u^ a u a is the four-acceleration of the fluid, 9 
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and u 2 = -oJ^uj 
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In a comoving reference frame the four-velocity of the fluid can be written as 
u M = [1,0,0,0]. The divergence of the four-acceleration will then vanish, and the 
Raychaudhuri equation can be simplified to 

q2 

= 6 + — + 2(a 2 - uo 2 ) + R 00 . (21) 
o 

Analogously to the Hubble parameter in FRW models, we can define an effective local 
Hubble parameter from the expansion rate of the fluid: 

H = . (22) 

This allow us to define an effective deceleration parameter in terms of H and H: 

q = = -1 - W/9 2 . (23) 

Using the relation R 00 = f (p + 3p), Eq. (|2*T|) can be written 

9 2 q = 6(a 2 -u 2 ) + ^K(p + 3p). (24) 

While Eq. (I24j) is valid for any metric, let us now specialize to the case of an LTB model. 
In it was claimed that the vorticity scalar is non-vanishing for the models considered 
in the paper, which may permit negative values of q. We find, however, that the vorticity 
tensor (|2()j) vanishes for a comoving fluid in the metric This incorporates also the 
models considered in For this metric the shear can be written as 

a 2 = ±(H r -H ± ) 2 . (25) 



Finally, assuming the fluid to be pressure-less, Eq. (j24j) reduces to 

9 2 q = 2(H r -H ± ) 2 + ^K P . (26) 

This equation shows that the dust dominated LTB-models have q > 0, meaning that 
they cannot have accelerated expansion. 

Although we have now shown that the LTB models cannot have accelerated 
expansion, this doesn't necessarily mean that they cannot explain the data used to 
infer this expansion. The first data used to conclude that the expansion seems to be 
accelerating were the measurements of the luminosities of supernovae of type la (SNIa) 
|17[ ITH] . Since then newer data have appeared that strengthen this claim even further 
|19| 12*0*1 |2*T] . However, accelerated expansion follows logically from these data only if 
one assumes that the universe is homogeneous. The added freedom of having a position 
dependent expansion in LTB models allows one to explain the data without the need for 
the expansion to accelerate locally. The explanation would then be that the expansion 
rate is highest at r = and decreases with distance from the center, since the oldest 
supernovae are also farthest away. In this case a decelerating universe would seem to 
be accelerating. Celerier shows in j22J how such inhomogeneities can mimic the effects 
of dark energy in supernova observations, when interpreted within the framework of 
FRW models. More specifically, Iguchi et al. |2Sj show that models with a(r) > can 
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reproduce the luminosity distance- redshift relationship of a ACDM model up to z ~ 1, 
but not for higher redshifts. However, when compared to actual SNIa observations, it 
is easy to find models that fit the data even better than the best-fit ACDM model [24J. 
In fact, as we show in [21], these models might also be compatible with the CMB power 
spectrum, even though they do not contain any form of dark energy. 

Finally, we wish to stress the point that in this work we have shown that the LTB 
models cannot have an expansion that is locally accelerated. This does not exclude 
the possibility that there can be a so-called volume averaged acceleration, where a scale 
factor defined via the physical volume of a comoving region is found to have a positive 
double time derivative. Recent papers discussing the possibilities of this effect explaining 
the apparent accelerated expansion are Refs. [23 123 12H1 123 EOl EH E21 IS3 EH EHJ 
At the present time there does not appear to exist a consensus among cosmologists on 
this issue. 

4. Conclusion 

Recently it was claimed by Kolb et al. pQ that inhomogeneous perturbations of FRW 
universe models can result in accelerated cosmic expansion, and thereby eliminate the 
need to postulate the existence of the mysterious dark energy. However, several papers 
have appeared since criticizing this work heavily. These attribute the effects claimed by 
Kolb et al. to incorrect handling of second order terms in the perturbative expansion. 
When dealt with correctly, one would then assume that these terms do not lead to 
accelerated expansion. However, this leaves open the possibility that the full, non- 
perturbative solutions of the Einstein equations for inhomogeneous models might exhibit 
accelerated expansion. 

In this work we investigated a special class of inhomogeneous models which can 
be solved analytically, the so-called Lemaitre-Tolman-Bondi (LTB) models. These are 
the spherically symmetric, inhomogeneous models. We found a general solution to a 
dust-dominated LTB model, and showed that there exist solutions that appear to have 
accelerated expansion. However, these solutions turned out to be unphysical. Finally, 
it was shown that any physically realistic solution of Einstein's equations must have 
decelerated expansion if it contains only matter. 
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